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ABSTRACT 

Using the local helioseismic technique of ring diagram we analyze the frequen- 
cies of high-degree f- and ^modes derived from both velocity and continuum in- 
tensity data observed by MDI. Fitting the spectra with asymmetric peak profiles, 
we find that the asymmetry associated with velocity line profiles is negative for all 
frequency ranges agreeing with previous observations while the asymmetry of the 
intensity profiles shows a complex and frequency dependent behavior. We also 
observe systematic frequency differences between intensity and velocity spectra 
at the high end of the frequency range, mostly above 4 mHz. We infer that this 
difference arises from the fitting of the intensity rather than the velocity spectra. 
We also show that the frequency differences between intensity and velocity do 
not vary significantly from the disk center to the limb when the spectra are fitted 
with the asymmetric profile and conclude that only a part of the background is 
correlated with the intensity oscillations. 

Subject headings: Sun: helioseismology -Sun: oscillations - Sun: interior 



1. Introduction 



Different helioseismic instruments, both from ground and space, observe the Sun in dif- 
ferent observables. Due to the different techniques used by these instruments, it is possible 
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to measure the solar oscillations simultaneously either as variations in the photospheric ve- 
locity or as intensity fluctuations. It is therefore important to confirm that the oscillation 
mode parameters measured from both the intensity and velocity agree with each other to a 
high degree of precision. However, the initial measurement of low degree p-mode frequen- 
cies from Doppler velocity (V^) and continuum intensity (Ic) observations from Michelson 
Doppler Imager (MDI) instrument on board Solar and Heliospheric Observatory (SOHO) 
showed systematic differences. A comparison of 108-day power spectra between V and Ic 
showed a weighted rn ean difference of —0.1 /xHz for i = 0, and —0.16 /iHz for 1=1 modes 



( iToutain et al. 1119971 ). Since the apparent frequency shift between an oscillation observed in 



velocity and intensity cannot be a property of the mode, it must arise from systematic errors 
while calculating the frequencies from the observed power spectra. Hence it was argued that 
the source of the systematic difference could be due to the opposite a symmetry effect seen 



betwe en the velocity and intensity power spectra (jPuvall et al. Ill993l ). lAppourchaux et al 



time 



(Il998f) also presente d a sim ilar evidence using VIRGO and SOI/MDI data. Around the same 



Harvey et al. I (119981 ) reported that the intermediate degree modes observed in V and 



total spectral intensity also show different central frequencies and observed that the apparent 
differences could be as large as 50 /iHz c lose to the acoustic cut-off frequency. However, the 



analysis of iNigam &: Kosovichev I (1l998l ) , using an asymmetric line profile-fitting formula, 
illustrated that the frequency difference between V and Ic in the intermediate degree range 
is much smaller compared to that obtained by fi tting a symmetric Lor entzian profile. Using 
the same asymmetric line profile-fitting formula, iToutain et al. I (119981 ) re-analyzed the data 
from MDI and showed that the frequenc y differences between V and Ic are considerably re- 
duced. iGavryusev fc Gavryuseva I (119991 ) have also analyzed data from different instruments 
and have argued that the reported frequency shift is merely an artifact of the reduction 
technique. 

Renewed interest in the topic began when local helioseismic te chniques were developed 



to stu dy the properties of high-degree modes in localized regions. iBasu. Antia. fc Bogart 



( 120011 ) analyzed azimuthally averaged (2-d) power spectra and inferred that the eigenfrequen- 
cies obtained using the asymmetric peak profiles agree reasonably well with each other while 
the use of symmetric profiles gives significant differences between frequencies computed using 
continuum intensity and velocity or continuum intensity and line-dept h spectra. In orde r 



to gai n further information for high-degree and high-frequency modes, iJain. Hill, fc Toner 



( I2OO3I ) analyzed the high-resolution GONG+ data. These authors also compared the az- 
imuthally averaged power spectra of a region centered on the equator and reported that the 
frequency dependence of the frequency shift between V and I is negligible below the acoustic 
cutoff frequen cy around 5.3 r nHz a nd substantial above the cutoff. These results supported 



the finding of iHarvey et al. I (119981 ). However, the conclusion is based on the visual compar- 
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ison of the peak frequency of the power spectra and may not necessarily be a true measure 
of the shift due to the reversal of the asymmetry and different background power between 
V and Ic spectra. 

It is now well established that line asymmetry of the solar power spectra alters the 
peak frequencies tha t are o btained under the ass umption that the lines are symmetric 



[e.g. lAntia fc Basu I (119991 ): iBasu fc Antia I (120001 ). However, the cause of the opposite 
asymmetry between the velocity and intensity spectra still remains inconclusive. The cur- 
rent understanding is that the reversal in the sign of asymmetry between the V and Ic 
spectra is due to the influence of the solar background noise that is correlated with the 



source ([Roxburgh fc Vorontsov 1 119971 : iNigam et al. I Il998l : ISeverino et al. I l200ll ) and the 



level depends on the characteristic granulat ion. On the other hand, the numerical simula- 
tion (iGeorgobiani. Stein, fc Nordlund II2003I ) indicates that the reversal is produced by the 
radiative transfer effects. Since the physics of the correlated noise is not yet fully understood 
and the spatial leakage signature for V and I is different due to their center-to-limb varia- 
tions, our objective is to examine the frequency dependence of the observed asymmetry and 
differences in eigenfrequencies between velocity and intensity observations as a function of 
the radial distance from the disk center to the limb. A prelimi nary investigation of a similar 
nature using azimuthally averaged power spectra is reported in lTripathy et al. I (120061 ). How- 
ever the present analysis differs from all earlier ones since here we use the three-dimensional 
(3-d) power spectrum, which is associated with flow fields, while the azimuthally averaged 
spectrum has no flow fields associated with it. 

The rest of the paper is organized as follows: We briefiy describe the data and analysis 
technique in Section 2, while the results are described in Section 3. Finally, we summarize 
the main conclusions in Section 4. 



2. Data and the Technique 



We use the Dopplergrams and continuum intensity images obtained by the MDI instru- 
ment during the period of 1997 May 19 - 21 when solar activity was near minimum. We 
have chosen 4 regions centered at heliographic longitudes of 0°, 15°, 30°, and 45°; all centered 
at the equator. The spatial extent of each of the localized region covers 256 x 256 pixels in 
heliographic longitude and latitude and gives a resolution of 0.01616 Mm"^. Each region is 
tracked for 4096 mi nutes , whic h gives a frequency resolution of 4.07 /iHz. The standard ring 



diagram technique (IHill Ill988l ) was used to obtain the power as a function of {k^, k. 



To extract the frequencies and other mode parameters, the three-dimensional power 
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spectrum is fitted with a model with asymmetric peak profiles of the form 



P{kxi ky, u) 



=B2 



+ 



k^ 



+ 



expjAo + {k- ko)A^ + A^i^)'' + Ag^)^, 



where 



X 



y,Vy 



wq + wi{k - ko) 
= S^ + {1 + Sx)\ 



(1) 

(2) 
(3) 



Here P is the oscillation power for a mode with a temporal frequency i' and the total 
wavenumber k = ^Jk'^ + k'^. The peak power in the mode is represented as exp(y4o) while 
the coefficients A\ to ^3 account for the variation in power with k and along the ring. The 
terms involving B\ and B2 define the background power and is the mode width. The 
term ck^ gives the mean frequency, and this form is chosen as it gives satisfactory fits to the 
mean frequency over the whole fitting interval. The parameters S c ontrols the asymmetry 



i n the peaks, and the form of asymmetry is the same as that used by iNigam fc Kosovichev 



( 1l998l ) i.e. the parameter is positive for positive asymmetry and negative for negative asym- 
metry. By setting = we can also fit symmetric Lorentzian profiles. A more detailed 
descri ption of the parameters and fitting procedure is described in iBasu. Antia. fc Tripathy 
(I1999I ) and iBasu fc Antia I (Il999l ). The 13 parameters (Aq, Ai, ^2, A3, B^-, c, p, S", U^-, 
[/,,, wn, w\] are determined by fittiri g the spectra using the maximum likelihood approach 
( lAnderson. Duvall. fc Jefferies 1119901 ). In this work we express k in units of so that k is 
identified with the degree I of the spherical harmonic of the corresponding global mode. 



3. Results and Discussions 

Figure [1] shows a characteristic i — v diagram obtained from the fits to the power spectra 
at the disk center using asymmetric profile corresponding to Eq. (1). In general we obtain 
about 750 and 600 modes from the velocity and intensity spectra respectively. Intensity 
spectra normally yields less number of modes due to the inherently lower signal-to-noise 
ratio. The number of fitted modes also decreases for regions that are away from the disk 
center. Comparing the number of modes between symmetric and asymmetric fits we find 
that the symmetric fit yields more modes. This implies that the inclusion of the asymmetry 
in the fitting formula reduces the number of modes that are successfully fitted. This could 
be due to some cross-correlation between S and other parameters of the model, particularly. 
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the background. Although, we fit the spectra up to radial order of n = 6, we restrict the 
analysis up to n = 4 modes since higher order radial modes have large estimated fitting 
errors. 



3.1. Asymmetry of Peak Profiles 

Physically, asymmetry is a result of an interaction between an outward-directed wave 
from the source an d a corresponding in ward-directed wave that passes through the region of 



wave propagation (IDuvall et al. 1119931 ). To illustrate the peaks in the power spectra and to 
visualize the asymmetry associated with different observables, we show an example of the 
azimuthally averaged power for t = 675 corresponding to the disk center spectra for V and 
Ic in Figure [2l It is apparent that the line profiles are asymmetric. At low frequency, the 
asymmetry agrees with the known results; the velocity spectrum has negative asymmetry 
i.e. more power on the lower frequency side of the peak while the continuum intensity 
spectrum show positive asymmetry i.e. more power on the higher frequency side of the 
peak. This reversal of asy mmetry is believed to arise from the addition of correlated noise 



to the amplitude spectra (INigam et al. Ill998l ). In addition, we notice that the asymmetry 



associated with velocity spectra appears to be reduced or even reversed at higher frequencies. 
Visual inspection of the azimuthally averaged power spectra show a comprehensible frequency 
shift of the central frequency between V and Ic above the cutoff frequency of 5.3 mHz. We 
also observe that the velocity power is higher near the disk center and decreases toward the 
limb in agreement with the predominately radial nature of the oscillatory velocity field. 

Figure [3] shows the asymmetry of the line profiles as a function of frequency obtained 
from fitting the three dimensional power spectra at disk center and 45° E longitude. As 
expected, this parameter is predominantly negative for all modes in the velocity spectra 
(upper panels) indicating that there is more power on the low frequency side of the peak. 
The asymmetry is minimum around 3 mHz while at low end of the frequency range the 
asymmetry of p2 and higher order radial modes is more likely to be zero within the estimated 
errors. Thus, in this frequency range the mode frequencies obtained from asymmetric and 
symmetric profiles should be nearly identical. Surprisingly, we d o not observe the asymm etry 



to be reversed at higher frequencies as noted above and also in iTripathy et al. I (120061 ). 



The asymmetry for the intensity spectra (lower panels of Figure [3]) appears to be more 
complex and in general these are larger than the velocity asymmetry. For / and ]?i-modes, 
this parameter is primarily positive (more power on the higher frequency end of the peak) 
and increases with frequency. For higher order p-modes, asymmetry is negative both at 
low- and high-end of the frequency range and is likely to be insignificant or zero at other 
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frequencies. Thus the asymmetry is different than what is obseryed for low-degree modes 



( iToutain et ah Ill998l ). intermediate degree modes (INigam et ah I Il998[). and high-degree 



modes obtained from the azimuthally averaged spectra (iTripathy et al. II2006I ). Ahhough 



most of these differences may be due to frequency, we do not rule out the possibility that the 
form of asymmetry used in the model are simplified and may not represent the real profiles 



adequ ately at all frequencies. The asymmetry may also weakly depend on I (IBasu fc Antia 



19991 ) and this dependence has been neglected in the fitting formula. Finally, we note that 



the asymmetry of both velocity and intensity spectra do not show any significant changes at 
different locations on the solar disk. 



3.2. Comparison of asymmetry between 2-d and 3-d fitting 

It is interesting to compare the asymmetries from 3-d spectra with those obtained from 
azimuthally averaged spectra. Figure H] illustrates the asymmetries obtained from both the 
fits. While the asymmetry observed in intensity spectra obtained from 2-d fits are entirely 
positive, the behavior is more complex in the 3-d spectra. On the other hand, asymmetry 
associated with 3-d velocity spectra is entirely negative while for 2-d spectra, it reverses 
sign at high end of the frequency range. Thus significant differences are seen between the 
asymmetry obtained from fitting the three-dimensional and azimuthally averaged power 
spectra. 

Although it is not easy to understand these differences due to the different functional 
form of the model profiles to fit different spectra (six parameters are fitted in azimuthal 
spectra as opposed to 13 in the three-dimensional spectra), we have made a preliminary 
attempt by constructing synthetic 3D spectra with S = —0.05. This spectra are then 
azimuthally averaged and fitted to calculate the asymmetry in peak profiles. Because of 
non-zero values of velocities {Ux, Uy) the process of azimuthal averaging also introduces 
some asymmetry and hence the value of S in azimuthally averaged spectra is not the same 
as that used in preparing the 3D spectra. The effect of varying the flow velocity {Ux & 
Uy) and background term, Bi is shown in Fig [5l The synthetic spectra were constructed 
with a uniform amplitude Aq = 50, width wq = 40 /iHz and included only the f-mode ridge. 
The background in synthetic spectra included only the first term with different values of Bi. 
The effect of velocity increases with k and hence the shift in asymmetry from input value 
also increases with k. Further, since Aq and Bi are constant, the peak to background ratio 
increases with k. It is clear that even when the flow velocities are increased beyond their 
normal values, the maximum change in S is about than 0.02 at the highest value of i and 
rather large velocity is required to make significant differences. Similar results are obtained 
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when spectra with different values of S were used. In all cases azimuthal averaging tends to 
increase the value of S, though the effect may be small. 



3.3. Source of the asymmetry 



Peak asymmetry is believed to have a contribution from the acoustic source loca- 
tion and a contributioii from the effects of correlated background noise from the source 
Roxburgh fc Vorontsov have argued that the line asymmetries occ ur due to the 



modes correlated with velocity background while to reach the same effect, iNigam et al. 



( 119981 ) have proposed a small amount of intensity-correlated background. The amount of 
coherent correlated background , as well as the need of this component to model both V and 
Ic , has also been debated (e.g. IStraus et al. I (120011 )). 



In order to understand the role of the background power in the reversal of the asymmetry 
between V and Ic , we plot the background power as a function of degree i of the mode in 
FigureO In the same plot, we also show the background power obtained from the symmetric 
fits (open symbols). It is clear that the background power obtained from symmetric and 
asymmetric fitting of the velocity spectra (upper panel) closely agree with each other. In 
contrast the intensity spectrum returns a significantly higher background power when fitted 
with the asymmetric line profile compared to the symmetric profile. This indicates that the 
line asymmetry is reversed in the intensity. Since granulation contrast decreases from the 
disk center to the limb, we also expect the background power to decrease near the limb. But 
no significant variation in power is noticed from the disk center to the limb which supports 
the idea that only a part of the background is correlated with the intensity oscillations 
(iNigam fc Kosovichev 1^9991 : ISarban fc Hiini2004[ ). 



3.4. Frequency Differences 

The frequency differences from the symmetric fits to both V and Ic spectra at four 
different longitudes are plotted in Figure [71 The frequencies obtained from the intensity 
spectrum are systematically higher than the frequencies from the velocity spectrum, as one 
would expect due to the reversal of the line asymmetry. For the disk center spectra, the 
differences, in general, are on the order of 10 fiB.z or less up to u < A mHz; the minimum 
difference of about 2-3 fiB.z is seen for p2 and higher radial order modes in the 5-min band. 
Since near the limb we primarily observe the horizontal velocity field which is dominated by 
granulation, it is expected that the frequency difference between intensity and velocity should 
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decrease near the limb due to a decrease in correlation with oscillations. But we find that 
the difference increases systematically from the disk center to the limb. This again supports 
the idea that only a part of the background is correlated with the intensity oscillations. 

The use of asymmetric profiles to fit the power spectra results in a better determina- 
tion of the mode frequencies and as a result the frequency differences between velocity and 
intensity is significantly diminished (Figure [8]). Similar to fits with symmetric profiles, the 
differences are mostly positive as the frequencies obtained from intensity fluctuations are 
higher than those from the velocity. However, the /-modes show some exception at the high 
end of the frequency where the peak frequencies observed in velocity are marginally higher 
(and clearly larger than the estimated errors) than those observed in intensity. The dif- 
ferences between higher order radial modes, mainly and p^, are slightly bigger but these 
modes also have a large uncertainty in the fltted frequencies as seen from the error estimates. 
Although subtle variations are evident between different longitudes, the frequency differences 
do not reveal any systematic trend between the disk center and 45° longitude. 

Although the asymmetric fltting formula reduces the frequency differences between the 
velocity and intensity spectra, there still remains some systematic differe nces of the order of 



5-10 mHz, which is much higher than the uncertainty in estimated errors. iBasu. Antia. fc Bogart 



( 120011 ) have reported a similar result based on the analysis of azimuthally averaged spectra 
and have argued that the assumed proflles does not adequately describe the velocity spec- 
tra where the signal-to-noise ratio is high. However, our result indicates that the form of 
asymmetry used for intensity spectra may not be a correct representation of the real proflles. 
To demonstrate this, we show the difference in eigenfrequencies between identical modes 
obtained from the power spectra at different longitudes with respect to the spectra at the 
disk center in Figure [91 The frequency difference in V shows little variation in longitude 
for frequencies less than 3.5 mHz. Variations between intensity spectra are higher than seen 
in velocity. Moreover, the shapes of the curve in Figure M for velocity are similar as are 
the shapes of the curves in intensity, leading to the frequency differences shown in Figure El 
Therefore it appears that the use of an asymmetric line proflle provides more accurate es- 
timates of the eigenfrequencies of solar oscillation for velocity rather than for intensity. In 
addition, fltting the intensity spectra is more challenging than fltting the velocity spectra 
due to low amplitudes and low signal-to-noise ratio. The number of modes fltted from the 
intensity spectra clearly supports this conjecture. 
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4. Conclusions 

Using the local helioseismic technique of ring diagram analysis we have studied the high- 
degree /- and p-mode frequency differences between velocity and continuum intensity data 
obtained from the MDl instrument during a quiet period. Since the spectra are known to be 
asymmetric, we fitted the three-dimensional spectra with an asymmetric model profile based 
on the known form of the asymmetry observed in velocity and intensity. The asymmetry 
obtained from fitting the velocity spectra agrees with the previous results while that of inten- 
sity shows a frequency dependent behavior. We further notice significant differences between 
the asymmetry obtained from fitting the three-dimensional and azimuthally averaged power 
spectra. 

Fitting the three-dimensional disk center power spectra with a symmetric Lorentzian 
profiles leads to frequency differences on the order of 10 /iHz or less up to u < 4 mHz 
between intensity and velocity. This difference is found to be increasing systematically from 
the disk center to the limb. The use of asymmetric profiles leads to frequency differences 
that are smaller than the differences resulting from the symmetric fits. However, systematic 
differences still remain at the high end of the frequency range, mostly above 4 mHz. We 
demonstrate that this difference arises from the fitting of the intensity rather than the velocity 
spectra. We speculate that the form of asymmetry used in the model for the intensity spectra 
is over-simplified and does not adequately represent the real profile at all frequency ranges. 

We also conclude that the frequency differences between velocity and intensity do not 
vary significantly from the disk center to the limb. However, variations between intensity 
spectra are higher than seen in velocity. This supports the idea that only a part of the 
background is correlat e d with the intensity oscillati o ns. In this context sev eral authors (e.g. 



Severino et al. I (120011 ): iBarban. Hill fc Kras I (120041 ): iToutain et al. I (120061 )) have advocated 



the importance of fitting the intensity-velocity cross spectrum along with a multicomponent 
model of the coherent background. This provides an impetus for future work using the 
high-resolution observations from Solar Dynamics Observatory. 



This study uses data from MDI/SOHO. SOHO is a mission of international cooperation 
between ESA and NASA. This work was supported by NASA grant NNG 5-11703 and NNG 
05HAL411. 
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Fig. 1. — The I — v diagram constructed from the asymmetric fits to disk center spectra for 
velocity (circles) and continuum intensity (pluses). 
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Fig. 2. — Normalized power for I = 675 corresponding to the azimuthally averaged and 
normalized velocity and continuum intensity disk center power spectra. 
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Fig. 3. — Asymmetry parameter S obtained from fits to velocity (top panel) and continuum 
intensity (bottom panel) spectra at two different locations. The left and right panels refer to 
the regions at disk center and at a longitude of 45° E, respectively. The symbols represent 
different radial orders and are explained in the right bottom panel. 
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Fig. 4. — Comparison of the asymmetry parameter obtained from fits to three-dimensional 
(filled symbols) and azimuthally averaged spectra (open symbols) at disk center. The top 
panel is for velocity and bottom panel is for continuum intensity. The errors are not shown 
for clarity. 
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Fig. 5. — The variation in the asymmetry parameter as a function of the flow velocity and 
background and is obtained from the fitting of azimuthally averaged synthetic spectra. The 
three values mentioned against each of the symbols represent t/^, Vy and B\ component, 
respectively. A value of 0.01 for and Uy corresponds to actual velocity of 43.72 m/s. 
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Fig. 6. — Comparison of the background power obtained from symmetric (filled symbols) and 
asymmetric (open symbols) fits to the disk center spectra. The upper panel is for velocity 
and bottom panel is for continuum intensity. 
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Fig. 7. — Frequency shifts between velocity and continuum intensity modes fitted using 
symmetric profile at four different longitudes. The locations are marked in each panel. 
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Fig. 8.— 



Same as Figure [7] but for modes obtained from the fit using asymmetric profiles. 




Fig. 9. — Frequency shifts between velocity (top panel) and continuum intensity (bottom 
panel) modes at three different longitudes with respect to disk center obtained from the fit 
using asymmetric profiles. The locations are marked in the top panel. 



